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ABSTRACT 


The dependence of the charge exchange lifetimes on the mirror 
latitude for ions mirroring off the geomagnetic equator has been 
re-computed using the improved hydrogen distribution models which have 
become available since the earlier calculation by Liemohn. The 
Chamberlain model, with the input parameters determined by recent satellite 
observations, has been used to define the spatial distribution of the 
neutral hydrogen environment through which the ring current ions traverse. 
The resultant dependence of the charge exchange lifetime, t, on mirror 
latitude, is best fit by the approximation =* Tg cos^'^Aj^^, where 
Tg is the charge exchange lifetime for the equatorial particles. This 
is a significant change from the cos^X^^^ approximation of the previous 
results. 


1 / 


The relative Importance of the charge exchange mechanism 1n the 
recovery of the geomagnetic storm time ring current process has been 
considered in numerous work since the early suggestion by Dessler and 
Parker (1959) that this mechanism could be significant. The lifetimes 
calculated by Liemohn (1961) for the tens of keV ring current protons 
had been the primary theoretical source for comparing observational 
results with the charge exchange mechanism ( Frank , 1967j Swisher and 
Frank , 1968; Russell and Thorne . 1970). As Smith et al. (1976) have 
pointed out, there now exists both more precise cross sections for 
charge exchange of protons and atomic hydrogen ( McClure . 1966), and a 
much improved hydrogen distribution model ( Opik and Singer . 1961; 
Chamberlain , 1963), which has been used by Proiss (1973) in computing 
the proton lifetimes. The dependence of the charge exchange lifetimes 
on the mirror latitude for ions mirroring off the geomagnetic equator 
has not been re-examined in light of the improved hydrogen distribution 
models, and even the most recent articles ( Tinsley . 1976; Lyons and Evans , 
1976) ^n the literature still employ the approximation given by Liemohn 
(1961) under the assumption of the neutral hydrogen density model of 
Johnson and Fish (1960). 

In this letter we use the Chamberlain (1963) model ( fipik and Singer . 
1961) to provide the spatial distribution of the neutral hydrogen 
environment through which the ions traverse. This model appears to be 
in good agreement with recent Apollo 16 Lyman alpha imagery of the 
hydrogen geocorona ( Carruthers , 1976) and Lyman alpha airglow measure- 
ments aboard several satellites ( Meier and Hanae. 1970: Bertaux. and Blamont , 
1973; Meier and Mange . 1973). The model assumes a spherically syirmetric 


hydrogen density and is characterized by three parameters; the tempera- 
ture and density at the exobase, and tne satellite critical altitude. 
Critical satellite altitude, is defined to be that altitude below 
which there exist closed orbit satellite particles, and above which the 
creation of satellite orbits is negligible so that the only orbits pre- 
sent are those that have perigees between the exobase and the critical 
satellite altitude. In order to determine the functional relationship of 
the charge exchange lifetime of ions mirroring off the geomagnetic 
equator, x , to the lifetime of equatorially mirroring ions, t„, we 
performed a full numerical integration for the mirroring particle motion, 
and fitted the computed ratio, T^j^/Xg, to a function of the form cos'^Aj^, 
where is the mirror latitude. The dipole geomagnetic field is assumed 
in our calculations. 

The mean lifetime, x^, for protons or any other species of ions con- 
fined to the equatorial plane for charge exchange decay with atomic hydro- 
gen is given. by 

T , 1 

e nlr^jov (1) 

where n{r^) is the hydrogen density in the equatorial plane, v is the 
ion velocity and a is the charge exchange cross section for that ion 
species. However, an ion oscillating between its mirror points, M and 
M', Figure 1, encounters a variable instantaneous density, n(r), of 
neutral hydrogen atoms ( Chamberlain . 1963} and, alsn, the time it needs 
to traverse an element of length along its guiding field line varies with 
the changing magnetic latitude. Since the ions make a large number of 
oscillations before charge exchanging with the neutral hydrogen, we can 
effectively average the instantaneous hydrogen density over a mirror path 
to get the time averaged density, n, that the particle encounters. This 
is expressed as 
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The above integration is performed along the magnetic field line 


r = rQ cos'^X 


(3) 


and so, 


ds = r^ cosXyj 4-3cos^X dX 


(4) 


The velocity parallel to the field line, V|| (r), is expressed 
for the dipole field ( Roederer , 1970) as: 


v,| (r) = V \j ^ (4-3cos^X)^ s in^ 


cos^X 


(5) 


where is the equatorial pitch angle and is related to the mirror 
latitude as 


sin^ 


cos®X_ 
m 

V4-3cos2x^ 


( 6 ) 


The mean lifetime, t^, of the mirroring particles then is given by 
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T = 

m 


n av 


( 7 ) 


We note that the denominator of equation 
the bounce period and its evaluation provides 
check on the numerical integration techniques 
The denominator integral to be evaluated thus 


(2) is one quarter of 
us with a consistency 
to be discussed later, 
is simply 


: ni P u 

) s f cosA (4~3cos X) dX 

° •' I _ sin^ an (4-3cos^X)^ 

. cos®X 


(8) 


This integral was first evaluated by Hamlin (1961) an,d it 
has a singularity when X=Xp^. The present calculation includes ^ 
derivation of this integral only as a check on the evaluation of the 
numerator in equation (2). If we introduce a new variable propor- 
tional to the magnetic field 

X = V±;3co A stn^ (gj 

cos°X 


equation (8) reduces to a form which can be implicitly represented as 



( 10 ) 
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Krylov (1962) has developed very powerful techniques for numerical 
integration using Orthogonal Polynomials. It can be shown ( Abramowitz 
and Stequn t 1965) that 


u 

j dy = w^ g(y^) + Rp 


( 11 ) 


where 


y^ = a + (b-a)x^ 


X, = 1 - 5,' 


positive root of Legendre polynomial P 2 p(x) 

w^ = 2 are the Gaussian weights of order 

pAn+l C(2n)!]^ (Or.) , 

^n “ 4n+T‘ [(?nTn^ ^ (n;, a < n < b . 


2n 


Thus, the integral is evaluated exactly whenever g(y) is a polynomial 
of degree (2n-l) or less. In our integration we use n=32. In Table 1 
we show the excellent agreement between the function f{aQ) calculated 
using the above method and the exact numerical calculation of Hamlin 
(1961). 
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Table 1 


Values of fla^) 

Mirror Latitude 
X m (degrees) 

0 

4- 

8 

12 

16 

20 

25 

30 

35 

40 

45 

50 

55 

60 


from Hamlin (1961) 

Haml i n 
Values 

0.7405 

0.7449 

0.7601 

0.7842 

0.8159 

0.8535 

0.9063 

0.9626 

1.020 

1.076 

1.129 

1.179 

1.224 

1.264 


Present Calculation 

Present 

Values 

0.740 
0.742 
0.757 
0 781 

0.i'12 

0.849 

0.902 

*0.956 

1.014 

1.060 

1.119 

1.169 

1.213 

1.252 
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The ratio t„/t can be computed for given mirror latitudes, X , 

Ml C Ml 

by calculating the two integrals in equation (2) to determine n, the 
time averaged neutral hydrogen density, and using equations (1) and 
(7). The calculated ratios are given in Table 2 for between 0° 
and 60°, and the corresponding equatorial pitch angle, a^, is also 
presented. The temperature at ;he exobase was taken to be T^, = 1000°K, 
the critical satellite particle altitude was taken to be R = 2.5R , 
and the L-value of the mirroring particles was L = 4. R_ is the geo- 
centric distance of the exobase taken at an altitude of 500 km. A 

least squares fit of the calculated ratio, to a function of the 

« 

form cos'^X^ determined a value of j = 3.4. This is in contrast to the 

previously used ( Liemohn , 1961) function cos X^^. 

The primary cause of the difference in the value of the exponent 

of the cosine function is the difference in the neutral hydrogen density 

distribution between the Johnson and Fish (1960) model which Liemohn 

(1961) used and the Chamberlain (1963) model which we have employed. 

In the Johnson and Fish (I960) model the hydrogen density tails off 

more rapidly as a function of radial distance than it does in che 

Chamberlain (1963) model (Figure 2). 

In Figure 3 we show the calculated ratio of lifetimes compared to 
3 4 

the fitted function cos X^^, and compared to the previously used ( Liemohn , 
1961) function cos°x^. As can be seen from Figure 3, the function cos’’* x^ 
is a good approximation to the calculated values. The difference between 
the two functional approximations is not as great for latitudes near 0°, 
corresponding to near 90° equatorial pitch angle, as it is at the larger 
mirror latitudes. The difference at 40° mirror latitude (21.6° equatorial 
pitch angle), however, represents nearly a factor of 2 increase in the 




Table 2 


Ratios of Charge Exchange Lifetimes for Mirroring Particles to Equatorial 
Particles for Various Mirror Latitudes 


Mirror 


Equatorial 

Latitude 

Ratio 

Pitch-angle 

Am (decrees) 

Ttn/-Te 

oa (decrees) 

0. 

1. 

90. 

5. 

0.988 

79.44 

10. 

0.954 

69.17 

15. 

0.898 

59.42 

20. 

0.823 

50.32 

25. 

0.731 

41.97 

30. 

0.628 

34.33 

35. 

0.519 

27.58 

40. 

0.411 

21.56 

45. 

0.309 

16.33 

50. 

0.218 

11.89 

55. 

0.142 

8.24 

60. 

0.084 

5.34 
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) 

charge exchange lUetimes for ions mirr’>r1ng at that latitude. 

As we stated earlier the hydrogen density in the Chamberlain 
model depends on three parameters; density, N_, at the exobase, tern- 

w 

perature T^, at the exobase and the critical satellite particle altitude. 

These parameters in turn vary with the geophysical conditions 

(solar activity, magnetic activity, season, etc,), and the MSIS model 

of Hedin (private communication, 1976) uses the measured geophysical 

conditions to determine N„ and T„. We have considered the effect on our 

c c 

functional approximation, varying these three para- 

meters. 

1) The ratio t„/t^ is independent of N. but the actual lifetime, 

T, does depend on the density at the exobase. 

2) The ratio of lifetimes does depend weakly on T_. The variation 

in for a temperature variation of 900°K to 1200°K is 

shown in Figure 3 at point (b) for = 30®, It should be I 

m I 

noted that, for the major geomagnetic storm periods during the 

early life of the Explorer 45 (S^-A) satellite ( Smith and Hoffman . 

1974), Hedin's MSIS model (private communication, 1976) determines 

T to fall in the range 970®K to 1050®K. In computing the ratios 

w 

presented in Table 2, T^ = 1000®K was used. 

3) The recent observations of the Apollo 16 Lyman alpha measurements 

are adequately described by a value of R__ = 2.5R^ particularly 
in the region near 3R^ ( Carruthers , 1976). However, for the 
purposes of examining the effect on we have considered 

a range of R__ values. The variation in t„/t_ for R.. between 

ill 6 
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2.25R and 2.75R is shown in Figure 3 at point (c) for « 40°, 

c c in 

The effects on Tj^/Tg of either a variation (900°K to 1200°K) 
or an variation (2.25R^ to 2.75R^) are quite small. We have computed 
T|^/Tg for different combinations of and R^^. in the stated ranges 
and the maximum variations in the ratio are shown in Figure 3 at points 
(a) for - 20° and = 50°, The ratio was again fitted to the 
cos'^X|^ function, with the points weighted inversely by the spread in 
T^/Xg at each and j « 3.5 was determined. In fitting the ratio 
using various extreme combinations we found, for T^ between 900°K and 
1200°K, for Rgg between 2.25Rg and 2.75Rg and for L-values between 2.5 

and 10, that j falls in the range j = 3.5 t 0.2, Recall, that for 

Tg = 1000°K, Rgg = 2,5Rg and L = 4, the functional approximation was 

Tm = T, cos^'^X„, 
me m 

'his new determination of the dependence of the charge exchange 
lifetimes ‘on the mirror latitude for ions mirroring off the geomagnetic 
equator requires that previous agreements or disagreements of the charge 
exchange mechanism with satellite observations need to be re-examined 
using at least the general relation 


T = t cos^’^ - X . 
me m 
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Figure 1 . 


Figure 2. 


Figure 3. 
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FIGURE CAPTIONS 


Cross section of earth dipole field in a magnetic meridian 
plane. 

Functional dependence of atomic hydrogen density on the 

radial distance from center of earth as given by the 

Johnson and Fish model (shown as solid line) and the 

Chamberlain model (shown as dashed line). Values of 

T » 1200°K and = 2.5R. are used. A unit density is 
c sc c 

assumed for each model at the exobase (500km). 

Ratios of lifetimes of mirroring particles to equatorial 

particles as a function of mirror latitude (and of equatorial 

pitch angle). Values shown by the x's correspond to input 

parameters of = 1000°K and R^^ = 2.5R^ in the Chamberlain 

3 4 

model. The approximate functional fit, cos ’ is shown 
by the solid line. For comparison, the function corresponding 

g 

to cos (Liemohn, 1961) is shown by the dashed line. 

Variations in the ratio for ranges of R^^ between 2.25R- and 

2.75R^ and T from 900°K to 1200°K are shown by points (a) 
c c 

(X^ = 20° and 50°). Variations in the ratios are shown for 
m 

fixed R = 2.5R and T ranging from 900°K to 1200°K at 

SC c c 

point (b), and for fixed = 1000°K and R^^ ranging from 

2.25R- to 2.75R- at point (c). 
c c 
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